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0^ ' The diffusion problem over a saddle is studied using a multi-dimensional Langevin equation. An 

0^ ' analytical solution is derived for a quadratic potential and the probability to pass over the barrier 

deduced. A very simple solution is given for the one dimension problem and a general scheme is 
' ■ shown for higher dimensions. 

> : 
o . 

■ I. INTRODUCTION 

^SJ I The Langevin equationjy has been applied to most fields of physics. It was solved several times for parabolic 
potential wells, see e.g. [^,|[r As analytical solutions can be derived for quadratic potentials only, it has really rather 
become widely used in numerical simulations. 

Our purpose is to establish an analytical expression for the diffusion over a potential barrier. In order to have a 
solvable problem, we assume that, around the saddle point, the potential can be approximated by quadratic functions. 
Some very simple expressions are obtained in one dimension. Because many processes obviously involve more than 
one coordinate, we extend our analysis to multi-dimensional cases. We thus derive an analytic expression for the 
distribution function of the Langevin equation, valid for multi-dimensional models, and then study the probability to 
over-pass the barrier. 

Our approach is only valid for classical diffusion satisfying the dissipation-fluctuation theorem. A solution for the 

■ one dimension Langevin equation in the ovcrdamped limit is derived in Ref. |^ for Levy flights but it cannot be 
I , . simply extended to multi-dimensional Langevin equations. 

O The interest of our approach will be shown in the case of heavy-ion fusion problems, for which Langevin equation 

^ r type simulations have been used by several groups The very small cross-section of such a mechanism makes 

^ numerical simulations very difficult, because very large statistics have to be computed. Our analytical expressions, 
^ . though using somewhat crude approximations, could be useful to extract some general trends. The problematics of 

' realistic calculations is out of the scope of this paper, where we only discuss general considerations concerning the 

. Langevin formalism. 

_ _ I 

II. THE LANGEVIN EQUATION 
A. Introduction 

To study the diffusion over a 1-D parabolic potential barrier, V{qi) = —raiUj\qi/2 with a given initial condition 
qio < and pio > 0, the Langevin equation reads, 

qi+Piqi-ujlqi^riit), (1) 

where ri(t) is a Gaussian stochastic force. As discussed in Appendix B, we rule out any anomalous diffusion process. 
The first two moments of this force are, 
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< ri(<) >= 



and 



nil 



(2) 



in agreement with the dissipation-fluctuation theorem. In the previous equations, T, mi and /3i denote the temper- 
ature, the mass and the reduced friction, respectively. All these parameters are assumed to be time, position and 
velocity independent, or at most very slowly varying, in the vicinity of the saddle. The symbol < > indicates an 
ensemble average. 

For any n-dimensional problem, one can generalize the previous approach, replacing the (qi, qi) variables by vectors 



+ f3Q-n^Q = R{t). 



(3) 



Such a canonical form of the problem results from two successive, very standard manipulations: 

- i) change the representation by transforming all tensors T, such as the friction tensor and the spring tensor, into a 
form M~^/^TM~-'^/^, where M is the usual mass tensor, naturally, 

- ii) change again the representation so that the spring tensor becomes diagonal. 

More precisely, assume some initial representation with a vector of degrees of freedom Z = {zi, ...Zn}, driven by a 
constant, symmetric mass tensor M, a constant, symmetric friction tensor 5, a constant, symmetric spring tensor S 
and a random vector force J^. The initial dynamical equation reads. 



Mz + gz -SZ = T{t) . 



This is equivalent to, 



M^f'^Z + M-^''^gM-^''^M^''^Z - M~^''^SM-^''^M^''^Z = M-^/^T{t) 



(4) 



(5) 



Let now U be that orthogonal matrix which lists the right eigenvectors of M ^^^SM as columns. Accordingly, 
= Un'^U-^, where fl'^ is diagonal. (Throughout this argument we rule out, naturally, those very 
exceptional cases where diagonalizations and/or inversions are singular.) Then the dynamical equation reads as well. 



u-^M^/^z + u^Hr^^^GM^^^^u U'^M^/'^z - u^Hi^^^z = u^Hr^^^Tit) . 



(6) 



With the definitions, Q = U-^M^/'^Z, f3 = U-^M-^^^QM-^^^U and R{t) = M-^/^T{t), this is nothing but the 
canonical form, Eq.(||). 

For 2-D for instance, we may obtain the first final degree of freedom as a "valley" direction to the saddle and the 
other final degree as a "confining" direction. 



Q 



Luf 

-iuj 



Pi P12 
P12 P2 



(7) 



and R(t) is a random force with two components. Such two components are correlated when /3 is non-diagonal but 
can be related, through a suitable matrix F, to a vector of independent random numbers. 



(8) 



with 



< vi{t)vi{t') >= 5{t-t'), < U2[t)v2{t') >= 5{t~t') and < i'i{t)v2{t) >=0. (9) 

The matrix F is real, but usually not symmetric. The fiuctuation-dissipation theorem, incidentally, which is easy to 
derive from the initial form, Eq.(||), states that FF-'" = 2Tj3, where T is the temperature and where the superscript T 
indicates transposition. All necessary details are found in Appendix A. 

It will be stressed again here that the matrices /?, and F take into account an overall multiplication of motion 
equations by the inverse of the (square root of the) mass tensor. It is easy to prove that this manipulation does 
not change the signs of the eigenvalues of the resulting matrices, and that such resulting matrices are usually non 
diagonal. In turn, the final spring tensor Vl^ can be made diagonal by the additional manipulation ii). 

The generalization to n dimensions is trivial. 
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B. Analytical solution 



Defining P — Q, it is easy to transform Eq.(|^) into a first-order differential equation, in a (2n)-D space, 



d 


'Q' 




■ 1 ■ 




Q' 




■ 0" 


di 


p 








p 


+ 


R 



(10) 



In tfie following we call "drift matrix" (DM) that block matrix T) = 



appearing in Eq.(p^. In the upper- left 



1 

-/3_ 

and upper-right corners of this 2n x 2n matrix, the symbols and 1 denote, respectively, the null and the unit n x n 
matrices, naturally. Whenever this DM can be diagonalized, and in the limit where the time and space derivatives 
of all physical parameters such as masses, drifts, frequencies, etc. can be neglected, the previous system can be 
transformed into. 



d 
dt 











= D 




_X2n _ 




_X2n _ 



(11) 



where D is the diagonal drift matrix. 



and 6 is the "rotation" matrix. 



D = e-^ve = 



ai 
.. 

a2n 



(12) 







- Xi - 




= 9 




Pi 






-Pn- 




-X2n- 



(13) 



This matrix 6 is the matrix of right (column) eigenvectors of the DM. Any normalization may be chosen for such 
eigenvectors. 

It can be stressed here that the eigenmodes Xi are linear combinations of both positions and momenta, hence the 
saddle dynamics should be visualized in phase space rather than coordinate space only. All such eigenmodes are 
expected to decay exponentially with time, except just one, corresponding to a resulting preferred direction along the 
valley, in phase space. 

In that same limit where the derivatives of the physical parameters (masses, etc.) can be neglected, such first order 
differential equations can be formally integrated into. 



■i{t) = XiC - Xio ^ dre-°'^ [aiiJ^i(T) -f ... -l-a„t/„(r)], i = l,2. 



, 2n. 



(14) 



where a is 



Here the 's are defined from the effects of both 9 ^ and T matrices (in matrix notation, a — 9 

a 2n X n matrix, and and T are n x n ones) and the i^^'s are uncorrelated white random numbers, extending Eqs.(|^) 
to n dimensions. Then, the Euler type variables {xi, ...,X2n) should have the same statistical properties as those of 



dr e 



[Q;iii^i(r) + ... + ainVnir)], 



dre 



[a(2n)lJ^l(T) + ... -I- a(2«)ni^n(T)]^ . (15) 



These are Gaussian random numbers whose correlations are easily evaluated from those of the t'i's, 



A,j{t) ^< Xi{t)xj{t) >^ / dre (°*+"^')^ (aaOfji + ... + ai„aj„) = 



I _ ^-t{ai+aj) " 



Ui + aj 



k=l 



(16) 



In a more compact matrix notation, 
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A{t) =< 



' xi{t) 

_X2n{t) 

Notice the occurrence of the matrix aa^ - 



Jo 



(17) 



qT\-1 



, which is trivially positive definite if 9 is real. 





2Tf3 

Using functional integral techniques, the full distribution function can then be easily evaluated (see Appendix B 
for details) and reads, 



W{xi, X2n,t; Xio, X(_2n)o) 



1 



Xl 
X2n 



(18) 



where the matrix elements of A{t) are those defined by Eq.(p^). This makes also a Gaussian distribution, naturally. 
If necessary, it is then easy to return to the original variables, (gi, qn',Pi, ■■■,Pn)- Since 



(19) 



where the second term of the r.h.s. is the average trajectory and the first one is the diffusion part, it is obvious that. 



"91 " 










Pi 


= e 




+ e 




-Pn- 




.X2n e'^'"*- 







- <qi[t)>- 




'qio' 


< qn{t) > 




qno 


<pl{t) > 




PlO 


-<Pn{t) >- 




-PnO- 



(20) 



and 



Ait) 



qi- < qi{t) > 

qn- < qnit) > 

Pl-<pl{t) > 



[qi- < qi{t) >, ..,9„- < qn{t) >,pi- < pi{t) >, ..,_p„- < _p„(t) >] 



.Pn- < Pn{t) > 

Jo 

Eventually, the full distribution function reads, 

1 1 





p 



(21) 



(22) 



W{qi, ..,Pn,t;qiQ, ..,p„o) 



(O \n rr^TTT^ ^'^P -obi- < qi{t) >, ■■,Pn- < Pn{t) >]A \t) 



qi- < qi{t) > 

Pn- < Pn{t) > 



(23) 



after renormalisation with the Jacobian. This result is well known, see e.g. 

C. Probability to pass over the saddle 

To evaluate the probability of passing over the barrier, we are interested in the "reduced" distribution obtained 
when all degrees of freedom but qi are integrated out. It is also necessarily a Gaussian distribution. 
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W{qi,t;qio, ...,g„o,pio, .-^Pno) 



■ exp 



(qi- < qi{t) >) 



21 



2<W 



The only remaining task is then to evaluate < qi (t) > and ct^j (i) . From Eqs (^,^2|) one gets 

< qi it) >= Xio e-^i* + ... + 0i(2„) X(2„)o e"-*, 



and 



< (i) - < (^11 ^1 e°^* + ... + ei(2«) a:^2„ e'^^"*)^ > 

2n 2n 2n 2n 



=1 j=l 

2n 2n 



1=1 j=l 
n 71 



l—l J — 1 v — 1 w—1 



2rEE' 



ai + aj 



(24) 

(25) 

(26) 
(27) 

(28) 



To go further and do physics, one needs the eigenvalues and the eigenvectors of the DM. This is not always feasible 
analytically for any dimension. However, the previous scheme can be applied to particular problems where the drift 
matrix is explicitly known. 

Let us first consider the 1-D and the 2-D cases, for which some general features will be shown. 

III. THE ONE DIMENSION PROBLEM 

The 1-D approach is interesting because of its simplicity. Intuitively, choosing as unique variable, the valley one 
and averaging all over the others should be enough in a first approximation. In this approach, the DM can be easily 
diagonalized and the diffusion energetical condition, time and probability can be easily calculated. 



A. Solution of the Langevin equation 

There are only one mass and one random force, mi and i?i, respectively. The latter is related to one random number 
vi, normalized according to Eq.(^. The corresponding matrix T boils down to one number only, which, according to 
the fluctuation-dissipation theorem, reads T = (2/3iT/toi)~^/^. The eigenvalues of the DM are a = {I3[ — /3i)/2 and 
h = -{(3[ + /3i)/2, with I3[ = {(3l + 4^^)^/^. Note that a > and 6 < 0. The matrix and its inverse read. 



1 



4 



(/3;-/?i)/2 (/31+/3i)/2 



The matrix a = 



is thus, 



a = (2T/3i/mi)i/2 



The eigencoordinates read, with pi = qi, 



' X' 




Y 





(/3i-/3()/2 1 
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.Pi. 



The Euler type variables, x and y, are then defined as. 



x{t) = Xe^''^ ~ Xo 
yit) = re-"* - 



and their statistical properties must be the same as those of 

(2r/3i/mi)i/2 



(/3i-/30/2 1 



dTe-''^i/i(T), / dTe-^^vi{T) 
Jo 



(29) 



(30) 



(31) 



(32) 



(33) 



5 



see Eqs.(h4]lq). The random number i^i{t) being Gaussian, x and y are also Gaussian random variables, with 



< ^2^t) ^(1 - < y^t) ^(1 - < x{t)y{t) >= 

ami 



bmi ' 



2Tf3i 



(a + b)mi 



I _ g-(a+b)t 



(34) 



see Eq.(p 



To evaluate the probability for passing over the barrier, one needs the following distribution function, necessarily a 
Gaussian in the present model. 



W{qi,t;qio,pio) = 



■ exp 



[gi- < qijt) >] 



According to the first row of 9, see Eq.(p9|), the valley coordinate is, in terms of the eigencoordinates. 



(35) 



(36) 



The first part corresponds to the diffusion and the second one to the average trajectory. It is trivial to obtain Xq and 
Yo from qiQ and pio according to Eq.(^, hence < qi(t) > . It is also trivial to obtain < ql{t) > from Eq. (|3^) . All 
told, elementary manipulations yield. 



< qiit) >= qioe-^''/' 



cosh (^^/3[?j + I sinh Q/?;t)] + 2 ^e-^^*/^ sinh Q/Jjt) , 



and 



<{t) =< ql{t) > - < qi{t) >2=- 



T 



miLof 



1 - e-ft* sinh^ + 1^ sinh ^ + 1 



These results are in agreement with the well-known calculation done for harmonically bound particles 
sign of the spring force should be changed. 



(37) 

(38) 
where the 



B. The critical initial kinetic energy 



Defining the critical energy as the kinetic energy, K = ^mip^Q, necessary to have half of the particles to pass over 

K,^[^l±^\ (39) 



the barrier, it is obvious that it corresponds to limj^+oo < qi{t) >= 0. From Eq.(|37|) it can be easily shown that 

where B — miw^^io/S is the barrier height. In the weak friction limit, /3i ~ 0, this condition becomes Kc = B which 
is a trivial result. In the overdamped limit, /3i 3> 2uji, it becomes, Kc = {(3i / loiY B . 

In the case of nuclei, using typical values, fiui = IMeV and f3i = 5.10^^s^^, the overdamped limit is usually a 
good approximation and a big kinetic energy is necessary to overpass even a very small barrier: Kc — lOB. 



C. Diffusion time 



As for many physical problems the diffusion process dynamically competes with some other processes, it is interesting 
to extract also the time necessary to reach the top of the potential barrier. When the previous condition, Eq.(^9|), is 
exactly fulfilled, then 

01+13', ^ 

<qi{t)>=qioe (40) 

The average trajectory exponentially converges to the top of the barrier with a typical time equal to 2/(/3i + /3^), 
which becomes 1/lui in the weak damping limit and l//3i in the overdamped one. 
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When the initial kinetic energy is higher than the critical one, Eq.(|39|), the average trajectory reaches the top of 
the potential barrier at Uop such as, 

»„,A,;,J.^f,/|_AV ,41, 



The previous equation becomes 



K 

coth(a;i ttop) = \ — in the weak damping case, (42) 



B 

lUJi 

2 ^^"7 \ ^ V s 



coth I i (3i ttop] — i %~\/ ^ 1 ) the overdamped one. (43) 



D. Passing probability 

In this model the probability at a given time that the particle has passed over the barrier is simply, 



P{t;qio,pio) 



■ erfc — 



W{qi,t; qiQ,pio)dqi 
<qi{t)>\ 



For large times (t 3> ^/P'l), 



<qi> 



V2CT,, ,/2i0i + P^P[) 



V2a,,{t) 



2uji 



(44) 
(45) 



— i-(-\/B — Vit^) in the undcrdampcd limit, 
Tpi 



IB uji jK 



in the overdamped one. 



(46) 
(47) 
(48) 



The passing probability is then known as a function of the initial kinetic energy and the temperature. It increases 
from to 1 around the critical value Kc when increasing the initial kinetic energy. The higher the temperature, the 
smoother is this increase. 



E. Heavy-ion fusion 

When two nuclei are colliding, the kinetic energy of the projectile should be higher than the contact energy derived 
by Bass [0 to observe some fusion events. This so-called extra-push energy is generally interpreted as an additional 
barrier due to nuclear forces that has to be overcome by the viscous nuclear matter . The critical energy derived 
above can then be seen as the extra-push energy here. 

To calculate the fusion probability of two cold colliding nuclei, the difficulty is then to evaluate the temperature. 
For particles that can reach the barrier top, at a distance R12, we will assume that all the remaining energy is totally 
dissipated. Therefore, with a level density aiev, we assume that aievT"^ = K — B, neglecting the collective kinetic 
energy. Then the fusion probability is. 



P{qio,Pio) = ^erfc 



aievB^ 

K -B 



1/4 




(49) 



In this formula, K is the available kinetic energy when the two nuclei are in contact, and B is the remaining barrier 
height: 
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K = Er 



E.r. 



(50) 



In this equation, Bsass is the Bass barrier ||ic| ] and Erot is the rotational energy, Erot = -E'cm^fmp/ ^12? where bimp is 
the impact parameter. 

Therefore, the fusion probabihty reads. 



P{E,m,b) ^^CTic 



1/4 



I 



F) — Bsass — B 



_2u^ 

f3i+l3[\ 



£^cm(l 



B 



(51) 



To compare this probabihty with other theoretical calculations, we need to choose parameters specific to our problem 
and deduce the potential and dissipation terms. The present paper being rather dedicated to a good qualitative 
understanding of the Langevin model, numerics with realistic parameters will be kept for a future paper. However, 
numerical simulations already available seem to indicate that at least two dimensions are necessary to really have a 
good understanding of the diffusion phenomena in this problem [p| . 



IV. THE TWO-DIMENSIONAL MODEL 
A. Passing probability 

Choosing qi as the valley variable, the associated distribution function, W{qi,t; gio, 920:PiO;P2o), is a Gaussian and 
the probability at a given time to pass over the saddle reads. 



IQ 




-P(i;<7io,g20,Pio,P2o) = / W{qi,t;qio,q2o,Pio,P2o)dqi (52) 

< qiit) >\ 



(53) 



The difficulty is then to find < qi{t) > and crq^{t), see Eqs.(p5|) and (p7|). When the two degrees of freedom are 
uncorrelated, namely when (3i2 = 0, the eigenvalues of the DM are, obviously. 



ai = (-/3i + ^Pf + 4cjj )/2 ( ai = (-ft + y /f3j + iu;j )/2 

a2 = (-/3i - y/W±My^ or J "2 = (-ft - y^F+MlZ? (^A) 

as = (-/32 + Vf3l - 4^| )/2 ] ag = (-ft + zy /-/?! + 4c.| )/2 ' ^ ' 

04 = (-ft - - 4'^l)/2 I "4 = (-ft - V-/3| + 4^l)/2 

if UJ2 < ft/2 or UJ2 > ft/2, respectively. Note that only ai is positive, the other eigenvalues are negative or have a 
negative real part. 

When P12 7^ 0, it is easy to show that there are always two real roots. The first one, oi, is positive and increasing as 
a function of (3i2. The other one, 04, is negative and decreasing as a function of the same. When the other two roots 
are real, it is also trivial to show that they remain negative. For large values of /3i2, which start reaching significant 
fractions of that maximum, (/3i/?2)^^^, which is acceptable for a semipositivity of friction, the other two roots may 
become complex conjugate. But the statement di{ai) < for z > 1, is still true. To illustrate the behavior of the 
eigenvalues, we consider the special case where f3i = f32 = I and Wi = W2 = to- The equal diagonal viscosities being 
taken as a unit, the polynomial equation whose roots are the eigenfrequencies of the problem reads, 

a'^ + 2a^ + ~ Lj'^ ^ (3^2^^ . (55) 

Fig.l shows, when /3i2 increases from to 1, the graphs of the four roots, or of their real parts when some of them 
become complex. Here to = 0.35, which corresponds to overdamping at the beginning, when /3i2 = 0. The merging 
of the two intermediate roots and their complexification are transparent. Once such roots have become complex 
conjugate, their common real part, however, remains negative. 
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0.5 
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-1.5 - 







0.25 



0.5 



0.75 



Pi 



FIG. 1. Real part of eigenvalues as functions of non diagonal friction. Overdamped regime, u = 0.35 < (3/2 — 0.5. Notice 
how two real roots fuse, then become complex conjugate, to generate one branch instead of two. 

Then Fig. 2 shows the case uj = 0.6, where the two intermediate roots are always complex, with a negative, common 
real part. 

0.5 

Re a 




-1.5 



1 1 1 1 1 1 1 1 1 








1 1 1 


1 1 1 1 ''"i*.. . 







0.2 



0.4 



0.6 



0.8 /3i2 1 



FIG. 2. Real parts of eigenvalues as functions of non diagonal friction. Underdamped regime, uj = 0.6 > (3/2 — 0.5. Since 
two eigenvalues are complex conjugate, three curves only appear. 



Hence, when t oo, only the "ai mode" survives. The leading terms in Eqs.(|25|) and (27) for long times are 



<qi{t)>=enXwe'''\ <jli{t) = 



9^1 e^°i*(aa^)i 
2ai 



(56) 



The result for t ^ oo reads. 



-f'(9io,920,Pio,P2o) ^erfc( 



All = (aa^)ii/(2ai). 



(57) 



It is similar to the result obtained in 1-D, provided we choose the proper coordinate. 
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Xi = i^jqi + aipi - (ii2 ^ p (-'^292 + aiP2)- (58) 

LU2+ ai+ [ha-i 

The condition to have half of the particles to overpass the barrier is then Xiq — 0. When the two degrees of 



freedom are uncorrelated {(3i2 — 0), one gets the same condition as in one dimension, Eq.(39). But in general, it is 
not possible to simply express it in terms of the initial kinetic energy because it also depends on the orientation of 
the initial velocity in the potential landscape. For the same reasons, the diffusion time is not easy to evaluate either. 
Note that in the previous graphs, (Figs. 1 & 2), ai is almost constant and could be approximated by the value given 
in Eqs.dl). 

Then we can evaluate An, 

All - -me-^fiz + P2{e-^fu + Wi2{e-^)i^{e-%^), (59) 

ai 

provided the diagonalisation of the drift matrix can be done. Note that only (6*^^)13 and {6^^)u occurs, which are 
coupling Xi with the velocity coordinates pi and p2- From Eq. (p8|), one gets, 

^2 ^2 2w| + 2a\ + /32ai 
[ljI + al+ /32ai)- 



All = Tai ( /3i - Phal :i^7Ji..^2 ) ■ (60) 



Again, when the two degrees of freedom are uncorrelated (/3i2 — 0) , one gets the same probability as in one dimension, 
Eqs.(|||6|). 



V. CONCLUSION 



In this paper, we showed a general scheme to solve multi-dimensional Langevin equations near a saddle point. In 
one dimension, the solution is very simple and the diffusion condition, time and probability can be easily expressed. 
This means that in the case of a simple one dimension model for heavy-ion fusion, one can analytically estimate the 
extra-push energy, the fusion time and its probability. 

For higher dimensions, the stochastic dynamics is again easily solved, but a difficulty remains, namely the explicit 
diagonalization of the drift matrix. In such conditions, a general analytic solution cannot be written, but the main 
features we found could be easily applied to very specific physical problems. In particular, we showed the difference 
between the dominant degree of freedom and the damped ones. The possible occurrence of complex eigenvalues 
corresponds to residual oscillations in subdominant degrees. 

In every case, the Gaussian solution to the problem contradicts the naive, intuitive expectation of a distribution with 
two peaks moving apart from the fusion barrier. In fact, the way over the barrier simply results from a competition 
between the drift of the center of the Gaussian and its spreading. 

For exotic noises leading to anomalous diffusion, such as Levy fiights, the Langevin equation can be analytically 
solved for a one dimension overdampcd motion where the Langevin equation reduces to a Smoluchowski one. In that 
case, the solution is given in Ref. 
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Appendix A 

In order to avoid possible confusions brought by the three successive changes of representation expressed by the 
matrices Af U and 0, we prove again the well-known fluctuation-dissipation theorem. Our precise purpose is to 
parameterize the random force R, see Eq.(^), in terms of dimensionless, independent random numbers Vi, normalized 
to unity, see Eq.(^). The starting point is a simplified form of Eq.(^), namely a situation where just the friction and 
the random force are present, 

MZ + gZ = T{t) . (61) 
This leads at once to a simplified form of Eq.(||), 

M^/'^Z + M-^l'^gM-^l'^M^I'^Z = M-^/'^T{t) . (62) 

It is convenient at this stage to diagonalize the real, symmetric, positive definite matrix B = M~^/'^QM~^/'^ , and 
obtain its representation in terms of its eigenvectors \i > and eigenvalues A^, 

n 

6 = ^|i>A. <i|. (63) 

i=l 

This decouples Eq.(|6^) into independent modes, relaxing separately towards thermal equilibrium for large times, 

< i\M^^^Z >= ~Xi < i\M'^/^Z > + < i\M-^/^T{t) >, (64) 

namely 

< i\M^'^Z > ->t^+oo /*dTe-^'(*-") < i|M-i/2_^(r) > . (65) 
Ja 

The left-hand side of this limit, Eq. (|65|) , can now be squared as a kinetic energy and then ensemble averaged to be 
equated with the Boltzmann energy T = (< i\M^^^Z >)2 . (We set the Boltzmann constant as a unit, fc = 1, as 
usual.) The right-hand side of Eq.(|6^), in turn, can be constrained by the following ansatz for ensemble averaging, 



< i|M-i/2jr(^) >< i|M-i/2jr(^/) > ^ s{t - r'), (66) 
where pi is an unknown normalization. The square of the right-hand side of Eq.(|65[) thus becomes, 



t — >OG 

rfre-2^-(*-) = ^, (67) 

hence the normalization result, pi = 2T\i. Let A be the diagonal matrix defined by the square roots (2rAi)^/^ of such 
normalizations pi. Let U be that orthogonal matrix which lists the eigenvectors \i > as columns. Let ly = {i^i, i^n} 
be a column vector of normalized, Gaussian, independent random numbers. Since our fiuctuation-dissipation theorem 
states that Af^i/^jp ^ J^^^i \i > {2T\i)^/'^ =UKv, the stochastic force R present in Eq.(|l) reads R — T v, with 
r = U~^UA. This explains why T is not expected to be symmetric. Eventually, the fluctuation-dissipation can be 
written in a more classical way, 

rr^ = 2r/3, (68) 



11 



where (3 =U -^^^QM -^^^U is the reduced friction matrix. 



Appendix B 



The solution of multi-dimensional Langevin equations is shown here for Markovian noises. 

Using the fact that the Euler type variables, (xi, . . . , X2n), have the same statistical properties as those of Eq.(p^ 
we find, 



W{xi, ...,X2n,t;Xio, ...,X2no) = < S[xi - Xi{t)]. . .S[x2n " X2n{t)] > 

dki dk2n 



2ti ■'■ 27r 



< exp i[ki, ..,k2n] 



xi - Xi{t) 
X2n - X2n{t) 



>, 



(69) 
(70) 



where the a;i(i)'s are shown in Eq.(14). The time integral will be discretized, with 5t — t/L, 



Xl{t) 
X2n{t) 



= hm ye-(^-i/2)^-D 



re St 



i/„((5t) 



(71) 



Here, D is the diagonal drift matrix, see Eq.(|l2|) and ly'jiSr) = J(^£_i-jg^ d^rvjij) are the results of random walks in the 
force space during the time interval 5t. Then the distribution function can be evaluated in its Fourier space, using 
the statistical properties of such random walks, 



W{xi, ...,X2n,t; Xio, ...,X(2„)o) 



dki dk2 



2tt 2tt 









Xl 


exp 


i [ki, . 


■, fen] 










_X2n_ 



p{ki, ...,fc2„,i). 



(72) 



with 



{ki, ...,k2n,t) = lim TT < exp -i[fci, .., fe,; 



-(e-l/2)STA 



> 



(73) 



We have used the fact that the random numbers are Markovian to say that the average of a product is the product 
of the averages of its factors. 

In the case of a Gaussian noise, the distribution function of the random numbers reads. 



1 



exp I [i^i, .., v„ 



(74) 



and, due to their Markovian properties, their integration over a time step 5t leads to 



p{vi{5T),..,Vn{5T)) 



Therefore, 



^(27r(5T)" 



5t 



exp I - — [vi{8t), .., i^„((5t) 



(75) 



(fci,...,fc2„,i) = lim TTexp -^[fci,..,fc2„]e-(^-i/2)^-^ aa^ e-^^-^/2)Sri 



fci 



(76) 



Reintroducing the time integral, one obtains, 



p(fci, fc2„,i) = exp I [fci, .., fc2„]A 



h 



2n 



(77) 



where the matrix A is defined by Eq.(17). Inverting the Fourier transform, upon p{ki, fc2„, t), one eventually gets. 
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W{xi, ...,X2n,t\X 



10, 



(2n)07 



(27r)" VdiFI 



exp --[xi, .., a;2,i]A ^ 



(78) 



For Levy flights, we will restrict ourselves to a one dimension problem. The noise is defined by its characteristic 
function p{k) in the Fourier space, 



p{k) = j dve-''"'p{v) =exphA|fc|^], 



(79) 



and is the source of an anomalous behavior characterized by a mean square displacement of the form < (Jq{t) >cx 2 AH 
[ p2[ . As the fluctuation-dissipation theorem is not satisfied any more, we keep a generalized diffusion coefficient, A. 
For the special case /x = 2, the noise is Gaussian. Similarly to the Gaussian case, the integration over a time step 5t 
leads to, 

p'{k) = J diyiST)e-''"'<-^^^ pii^idT)) =cxp[-A{STy-''\k\''], (80) 
after renormalization ||l3|. Therefore, the average value in the r.h.s. of Eq.(|7^) can be calculated, 

Reintroducing the time integral, one finally gets, 

-A / (fcie-"^ + fc2e-''^)''dr 

"'0 



p(fci, fc2,t) = exp 



(82) 



For ^ — 2, i.e. in the Gaussian case, the time integral can be evaluated analytically, but such is not the case for any 
value of /i. The only favorable situation occurs for a one dimension Smoluchowski equation where. 



p{ki , t) = exp (^- A /cf e^^^^dr^ , 
see Ref. [|. This is why we restrict our study to classical Gaussian noises. 



(83) 
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